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Nonfactorizable effects in B to charmonium decays 
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Abstract. Nonleptonic B to charmonium decays generally deviate from the factorization predictions. We 
study rescattering effects mediated by intermediate charmed mesons in this class of decay modes and, in 
particular, we consider B~ — > K~h c with h c the J PC — l" 1 cc meson, relating this mode to B~ — > K~XcO- 
We find B(B~ — > K~h c ) large enough to be measured at the B factories, hence this process could be used 
to study the poorly known h c . 
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1 Introduction 

Testing the Standard Model description of CP violation in 
the B sector requires a reduced theoretical uncertainty, a 
difficult task for nonleptonic decays since a general compu- 
tational scheme has still to be developed. The amplitude 
of two-body nonleptonic B decays is given by the matrix 
element of the effective Hamiltonian for B — ► M\Mi HH: 



G 



A(B -» MiM 2 ) = ~^Y. ^Ci{ji){MiM2\Oi{ji)\B) 



(1) 



Xi are CKM matrix elements, Ci(fi) Wilson coefficients 
evaluated at the scale fi and Oi four-quark operators. The 
naive factorization ansatz expresses their matrix elements 
as products of matrix elements of quark currents. 

Let us consider B~ — > K~M 5c , where M 5c belongs to 
the charmonium system. Neglecting the annihilation term 
(suppressed by V u b), the factorized amplitude is: 



A F (B- -> K-M- cc ) = -j=V cb V* a 2 (n) + £ OiQi) 

V 2 - i=3,5,7,9 

(K-\(sb) V - A \B-) x {M- cc \(cc) Vt a\0) (2) 



with ci2 — C2 + Ci/N c , di = Ci + Ci+\/N c . Eq.Q shows the 
drawbacks of the approach: scale and scheme dependence 
of the Ci(fi) are no more compensated by that of the matrix 
elements; besides, the product of these does not contain 
any strong phase. An improvement consists in the general- 
ized factorization, with ai(/j,) replaced by effective (process 
independent) parameters af ff to be fixed using experimen- 
tal data. Other methods, like QCD-improved factorization 
El, PQCD 0, SCET g|, QCD sum rules 0, can only be 
applied to selected classes of nonleptonic modes. 

Generalized factorization indicates nonfactorizable ef- 
fects in B — > K~J/ip. The experimental B(B — ► K~J/i/j) 
can be fitted with \a^\ — 0.2 — 0.4 depending on the form 



factor parameterizing (if ~|(s&)v-a|-B~); |a2 ff | = 0.38 ± 
0.05 is obtained using the form factor in This must be 
compared to 02 = 0.163(0.126) computed in the naive di- 
mensional regularization (or 't Hooft-Veltman) scheme . 
The difference between a 1 ^ and 0,2 signals nonfactorizable 
effects. However, the clearest evidence of deviation from 
factorization is the observation of B~ — > K~XcO, with Xco 
the lightest cc scalar meson. The data: 



B{B~ -> K - Xc0 ) 
B(B- K- Xc0 ) 



(6.0 



+2.1 
-1.8 



± 1.1) x 10" 



(2.7 ±0.7) x lo- 



ft 
(4) 



from BELLE and BABAR respectively, show that 
the experimental amplitude is non-zero, while the factor- 
ized amplitude vanishes since (Xco\(cc)v^a\0) = 0. Be- 
sides, the rate is comparable to B~ — > K~J/4> since, for 

example, gg^O^ = (0.60t°;?J ± 0.05 ± 0.08) Q. 

QCD-improved factorization does not reproduce the 
measured branching ratios for B~ — > K ~Xc0, K~ J /if), giv- 
ing either small results or producing infrared divergences, 
a signal of uncontrolled nonperturbative effects UJ. 

In |10| we investigated if the deviation from factoriza- 
tion in B — > cc decays may be ascribed to rescattering of 
intermediate charm mesons as in the diagrams in figQ We 
found that such effects could be large enough to produce 
the observed B(B~ — > K~Xco)- Other modes with vanish- 
ing factorized amplitude can test the rescattering picture, 
as B~ — > K~h c with h c the lowest lying J PC = 1 H cc 
state. h c is listed by the PDG among the particles to be 
confirmed, with m hc ~ 3526 MeV 11 . If B~ -> K~h c 
has a sizeable rate, this decay may be used to study h c . 

To improve the analysis in |10| we introduce an ef- 
fective lagrangian for the interactions of all the low-lying 
(*) 

i — 1 cc states to mesons, based on the spin sym- 
metry for the heavy quark in the infinite mass limit. This 
relates all the couplings to a single parameter. The same 
holds for the couplings of t = cc mesons to D K ,J . 
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Fig. 1. Typical rescattering diagrams contributing to B 
K~M C c, with Mcc belonging to the charmonium system. 



2 Calculation of rescattering diagrams 

The factorized amplitude Ap(B~ — > K~h c ) in J5| van- 
ishes since (h c \{cc)v^A\fy = due to conservation of par- 
ity and charge conjugation. However, the decay can pro- 
ceed by rescattering induced by the same (bc)(cs) effective 
Hamiltonian. We consider the process B~ — > X^Y^ — > 
K~h c . The lowest lying intermediate states X® c , are 
Dg~, Z)(*) rescattering by exchange of To com- 

pute the diagrams in fig^we need the weak vertices B — > 
£)(*)£)(*) and two strong vertices: the coupling of charmed 
mesons to kaon, and of the h c to mesons. All those 
vertices are related to few parameters when wlq — > oo. 

Interactions of mesons Hq with a single heavy quark 
Q can be described using the Heavy Quark Effective The- 
ory, exploiting the heavy quark spin and flavour symme- 
tries holding for rriQ — > oo. In this limit the heavy quark 
four velocity v coincides with that of the hadron and is 
conserved by strong interactions J5]- Because of the in- 
variance under rotations of the heavy quark spin sq , states 
differing only for Sq are degenerate and form a doublet. 
When the orbital angular momentum of the light degrees 
of freedom relative to Q is £ = 0, the states in the dou- 
blet have J p = (0~,1~), corresponding to (£>( s ), Dj*^), 
(-B( s ) , B^ ) . The doublet is represented by the matrix H a = 
P + [M^7 Jli -M„7 5 ] with P± = (1± ^)/2. is the vector 
state, M the pseudoscalar one (a is a light flavour index). 
M a , M* contain a factor \/m, with m the meson mass. 



Let us consider B 



D (*)- D {*)0^ for which factor _ 

ization empirically works |13| . The factorized amplitude 
is: 

{D^-D^\H W \B-) = ^V cb V c * sai 

(£>W°|(y - AY\B-) x (DW-|(V-A) M |0) (5) 

with cii — ci -\-C2/N c . In the heavy quark limit, the matrix 
elements in (0 can be written in terms of the Isgur-Wise 
function £, and a single leptonic constant F |12| . 

The '£)(*) K couplings, in the soft kaon limit, are 
related to a single parameter g through the effective la- 
grangian describing strong interactions of charmed mesons 
with the octet of the light pseudoscalar mesons ^1]: Cj = 
i g Tr[H bl ^AlH a ]. In £ /: A^ a = \ (tfd^ - ^) ba , 

H a = j Hl"/° and £ = e¥~ , with f^U = 131 MeV and 
M. a 3 x 3 matrix with the light pseudoscalar meson fields. 

We now consider the strong vertex involving the h c . 
For Q1Q2 mesons heavy quark flavour symmetry does not 



hold any more, but degeneracy is expected under rota- 
tions of the two heavy quark spins. This allows us to 
build up multiplets for each value of £. For £ — one 
has a J p — (CP, 1~) doublet. The corresponding matrix 
is [13: RiQiQ*) = P+lL^f, - L 75 ]P_ with = J/ip, 
L = r\ c for cc. For £ = 1, the multiplet is pWiSMA* = 

In the case of cc, xi = Xc2, Xi = Xci, Xo = Xco correspond 
to the spin triplet,_/ii = h c to the spin singlet |16| . 

For £ = 1 Q1Q2 states, the most general lagrangian 
describing the coupling to two heavy-light mesons Qiq a 
and q a Qi can be written as follows |17|: 



L\ = i—l r 



p(QiQ^H 2altl H la 



h.c. + (Qi <-> Q 2 



where H 2a — [M'al^ — M' a j 5 }P_ describes q a Qi mesons. 
C\ is invariant under independent rotations of the spin of 
the heavy quarks. g\ describes the interaction of heavy- 
light mesons with the three Xc states and with h c , relating 
the couplings in absolute value and in sign as well. This 
allows a proper analysis of the diagrams in fig.l. 

The interactions of £ = states with the 

heavy-light J p = (0 _ ,1 _ ) mesons proceed in P-wave and 
can be described by a lagrangian with a derivative: 



r 92 rp 

£2 = -rTr 



R(QiQ2)H 2a @ H la 



h.c. + (Qi <-> Q2) 



again invariant under heavy quark spin rotations. 

gi and g 2 can be estimated using vector meson domi- 
nance (VMD). One can consider (D(v')\cc\D(v)), assum- 
ing the dominance in the t-channel of the + cc state, and 
using the normalization £(1) = 1. One obtains goDx^a = 
2m D m Xc0 /f Xc0 , where f Xc0 is defined by (0|cc|Xco(g)) = 
fxao m xco- Tni s relation gives: 51 = -y/m^/(V3f Xc0 ). 
The same argument gives 52 = y / m- l f,/(2mD fif>), where /y, 
is defined by (0|c7 M c| J/ip(p, e)) = f^m^e^. To account for 
the off-shell effect of the exchanged particles, the virtuality 
of which can be large, we write: gi(t) = gio Fi(t), with gio 
the on-shell couplings and: F l {t) = (Af — m 2 D( , ) )/(Af — t). 
The parameters Ai represent a source of uncertainty. 

In the numerical analysis we exploit the heavy quark 
limit, putting f D * = f Ds and use f Ds = 240 MeV 0. 

For the Isgur-Wise function, we use £(y) = (2/(1 + y)) 2 , 
compatible with data from semileptonic B — > Z)W decays. 
As for g, CLEO Collaboration obtained: g = 0.59 ± 0.01 ± 
0.07 by measuring r(D*) and the D* branching fraction 
to Dtt PHI- This value should be compared to theoretical 
estimates ranging from g ~ 0.3 up to g ~ 0.77 JHI- Since 
the rescattering amplitudes always depend on g ■ Fi(t), we 
put Ai = A for all i and use g — 0.59, leaving to A the 
task of spanning the range of possible variation of g. 

For (71 and 32 we use the VMD relations, together with 
the QCD sum rule result / Xc0 = 510 ± 40 MeV [ID] and 
the experimental value /ju = 405 ± 14 MeV. 

We computed the imaginary part of the rescattering 
diagrams; the determination of the real part is more un- 
certain. Using a dispersive representation, we obtained for 
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B~ — > K~Xco{J/i)) that KeAi ~Im„4i |10| - Hence we ac- 
count for the real part of the amplitudes considering them 
as fractions of the imaginary part varying from to 100%. 
Such an uncertainty will affect the final result. 

We constrain A considering rescattering contributions 
to B~ -> K~J/tp, where A{B~ -> K~J/ip) = A fact + 
Aresc is bounded by experimental data. For A ~ 2.6 — 3 
GeV, the sum does not exceed the experimental bound. 
Moreover, we repeat the analysis in |E3] for B~ — > A~XcO: 
using now the relations stemming from the lagrangian C\. 
We get a branching ratio compatible with the experimen- 
tal result from BABAR if A is varied around 3.0 GeV. 

Using such constraints we analyze B~ — > K~h c . In 
figGlwe plot the branching ratio versus A, We find |17j: 



B(B~ 



K~h c ) = (2- 12) 



1CT 



(6) 



where we accounted for the uncertainty on the real part 
of the amplitudes and on the variation of A. This result 
suggests B(B~ — > K~h c ) large enough to be measured 
at the B-factories. Moreover, this mode represents a size- 
able fraction of the inclusive B~ — > X/i c decay, for which, 
considering the production of the cc pair in /i c in the color- 
octet state, B(B~ -> /i c X) = (13-34) x 10~ 4 is predicted 

EH- 

The main uncertainty affecting our results is due to 
cancellations between different amplitudes, which are of 
similar size. Another uncertainty is due to the neglect of 
contributions of higher states, even though a minor role 
can be presumed for higher resonances. 

Bearing such uncertainties in mind we can conclude 
that rescattering terms may contribute to the nonfactor- 
izable effects observed in B — > charmonium transitions. 



3 Conclusions 

The h c was observed in pp annihilation and in p — Li in- 
teractions I n B decays, one could access h c look- 
ing either at its hadronic modes: h c — > J/tpir , p°ir°, 
etc. , or at its radiative modes: h c — > ?y c 7, x c o7, etc. 
The channel h c — > 7^7 seems promising, as noticed by 
Suzuki who estimated: B(/i c -> r? c 7) ~ 0.50 ± 0.11 
A similar result: S(/i c — * ??c7) = 0.377 is obtained in 
|23| . These two predictions, together with the experimen- 
tal B(rj c — > KKtt) and our result JJjJ, allow us to predict: 
B(B- -> ii:-/i c -► AT"77 C 7 -► K-KKwy) = (4 - 26)10- 6 , 
within the reach of current experiments. 

As for rescattering effects in _B — > charmonium de- 
cays, we found that they can produce a branching ratio 
for B~ — > K~h c comparable with that of B~ — > A~x c o. 
The same holds for i? - — > A~V(3770) which, because 
of the smallness of ^(3770) , is predicted by factorization 
with a tiny branching ratio. The experimental measure 
A'-^(3770)) = (0.48±0.11±0.12) x 10^ 3 |H] is a 
further evidence of large nonfactorizable effects. In our ap- 
proach, using 3dd?/;(3770) = 14.94±0.86 obtained from the 
width of ^(3770), we get B(B~ -> K~iP(3770j) = (0.9 - 
4) x 10~ 4 . Similar conclusion applies to B~ — > K~\c2 
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Fig. 2. Branching ratio B(B~ — » K~h c ) versus yl. The lowest 
curve corresponds to Re^li = 0, the highest one to Re^li = 
Im^li. The dark region corresponds to the result @. 

with Xc2 the J PC = 2 ++ cc state, the factorized ampli- 
tude of which also vanishes. The observation of this pro- 
cess with branching ratio comparable to B(B~ — > A _ Xco)i 
B(B~ — * K~h c ) would support the rescattering picture. 
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